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^^stratospheric  behavior,  in  combining  theory  and  actual  measurements,  is  e power- 
ful diagnostic  tool, in  the-study  of  the  atmosphere^-  A model  is  very  useful  as 
an  aTdin  developing  both  the  general  approach  and  the  details  of  a field 
measurements  program,  and  as  an  important  tool  in  interpreting  the  experimental 
data  of  such  a program.  The  interplay  between  measurements  and  model  should 
produce  the  most  effective  approach  to  study  of  the  stratosphere. 

A few  examples  of  calculational  results  from  the  ASL  Numerical  Model  of  Atmo- 
spheric Radiation  (ANMAR),  Composition  and  Dynamics,  paralleling  the  conditions 
of  the  recent  set  of  STRATCOM  VI  experiments,  have  been  presented  toydemonstrate 
the  range  of  results  obtainable  and  the  detail  of  treatment  possible  using  the 
modeling  approach. 
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INTRODUCTION 


The  chemical  composition  of  the  stratosphere  is  an  ever-changing,  highly 
variable  result  of  the  complex  interplay  of  solar-induced  photodissocia- 
tion processes,  resultant  photochemical  and  thermochemical  reactions, 
and  a variety  of  transport  processes  on  space  scales  ranging  from  molec- 
ular to  global.  The  magnitude  or  relative  importance  of  any  process 
in  determining  composition  and  atmospheric  behavior  is  dependent  on 
altitude,  latitude,  time-of-day,  season-of-year , and  other  variables. 
Because  of  the  great  complexity  and  variability,  the  fundamental  require- 
ment for  gaining  understanding  of  the  stratosphere  is  measurement  of  a 
number  of  critical  atmospheric  parameters,  including  particle  densities 
of  key  trace  constituents,  solar  flux  and  thermodynamic  properties. 

These  measurements  in  themselves,  however,  do  not  provide  explicit 
information  of  the  processes  which  produced  the  results  observed.  They 
allow  prediction  of  the  state  or  behavior  of  the  stratosphere  only  in  a 
statistical  sense  if  large  numbers  of  measurements,  under  a variety  of 
time  and  space  conditions,  are  obtained. 

A supplemental  (and  necessary)  parallel  approach  is  the  development  of 
mathematical  simulations,  i.e.,  "models,"  of  stratospheric  composition 
and  behavior  based  on  the  physical  laws  which  apply.  Field  measurements 
and  theoretical  speculation  are  the  usual  modes  of  motivation  for  model 
development,  and  laboratory  measurements  provide  the  quantization.  A 
model  is  modified  as  required  by  the  accumulation  of  relevant  field  mea- 
surements; and  ultimately  the  degree-of-val idity  of  a model  is  determined 
by  comparison  of  calculated  results  with  actual  atmospheric  measurements. 
The  relative  validity  establishes  the  degree  of  confidence  that  can  be 
placed  on  the  diagnostic  or  predictive  capabilities  of  the  model. 

During  the  past  8 years,  the  Atmospheric  Sciences  Laboratory  (ASL)  at 
White  Sands  Missile  Range,  NM,  has  developed  a large,  complex  program 
of  balloon-based  stratospheric  measurements  in  cooperation  with  several 
other  laboratories.  One  of  the  central  themes  of  this  STRATCOM  (STRAT- 
ospheric  COMposition)  program  [1]  is  simultaneous  measurement  of  sets 
of  related  composition,  thermodynamic  and  radiative  parameters.  Such 
measurement  of  related  parameters  under  the  same  conditions  allows  direct, 
rather  than  statistical,  analysis  and  interpretation  of  relationships. 
Development,  testing  and  validation  of  computer-based  models  of  the 
atmosphere  can  also  be  more  direct,  and  confidence  in  their  capabilities 
is  enhanced  by  relating  to  actual  measurements. 

To  support  and  supplement  this  field  measurements  program,  a chemical- 
kinetic  model  of  the  stratosphere  [2],  incorporating  a parametric 
application  of  the  vertical  transport  processes,  has  recently  been  made 
operational  at  ASL.  The  model  is  used  as  one  of  the  tools  in  developing 
the  overall  measurements  program,  as  well  as  in  designing  individual 
experiments.  It  will  also  be  used  as  an  exploratory  tool  in  theoretical 
studies  to  extend  the  capability  of  interpreting  the  experimental  results 
obtained. 
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This  document  presents  examples  of  the  kinds  of  calculational  results 
obtainable  from  the  model. 


FIELD  MEASUREMENTS 


AND  ATMOSPHERIC  MODELS 


A recent  operation  in  the  STRATCOM  atmospheric  measurements  program, 
STRATCOM  VI  [3]  was  held  in  late  September  1975.  The  principal 
set  of  measurements  was  in  the  altitude  range  of  25-39  kilometers  during 
a 34-hour  period  on  24  and  25  September.  A supporting  set  of  IR  absorp- 
tion measurements  was  made  from  an  altitude  of  31  kilometers  at  sunset, 
26  September  1975,  using  a second  balloon.  Both  flights  were  from 
Holloman  Air  Force  Case,  NM,  32°  N latitude^ 


To  parallel  these  field  measurements,  the  calculations,  from  which  the 
following  sample  results  are  taken,  were  made  by  use  of  photolytic  data 
derived  from  solar  conditions  for  the  latitude,  dates,  and  altitudes 
given  above.  The  photodissociation  coefficients  were  calculated  by 
0.  L.  Collins,  and  are  presented  and  discussed  in  [4]. 


THE  ATMOSPHERIC  CHEMICAL-KINETICS  MODEL 

The  structure  and  calculational  method  used  in  the  computer  simulation 
of  stratospheric  composition  is  discussed  in  greater  detail  in  [2].  A 
brief  summary  is  given  here. 

The  computer  model  uses  the  Gear  method  of  solution  of  sets  of  stiff 
differential  equations  [5]  to  solve  the  coupled  continuity  equations 
describing  the  time  dependence  of  the  particle  densities  of  30  chemically 
reactive  atmospheric  species.  (Oxygen  and  nitrogen  molecules  are  included 
with  suitable  steady-state  densities.)  The  atoms,  molecules,  and 
radicals  treated  are  given  in  Table  1. 

! TABLE  1 

ATMOSPHERIC  CHEMICAL  SPECIES  CONSIDERED  IN 
[ COMPUTATIONAL  MODEL 

» 

0 0(iD)  0(lS)  O2  02(^4)  02(1e)  O3 

i N N2  NO  NO2  N2O  NO3  N2O5 

i 

H H2  OH  H2O  HO2  H2O2 


CO  CO2  CH2  CH3  CH4  CHO  HCHO 


A total  of  34  photodissociative  processes  and  115  chemical  reactions, 
coupling  the  densities  of  the  species,  are  considered.  These  are 
tabulated  in  [2]  and  partially  sketched  here  in  Figs.  1,  2,  and  3,  which 
show  families  of  atmospheric  constituents  and  the  processes  of  conversion 
from  one  species  to  another.  Figure  1 is  for  part  of  the  nitrogen/oxygen 
fami ly : 


0,  O2,  O3,  NO,  NO2,  N2O. 

In  the  diagram,  a chemical  species  at  the  tail  of  an  arrow  is  converted 
to  the  one  at  the  head  by  interaction  with  any  species,  or  a photon  (hv) 
listed  on  the  arrow  shaft. 

Figure  2 is  a similar  representation  of  some  of  the  reactions  of  the 
oxygen/hydrogen  family: 


OH,  H2O,  HO2,  H2O2 

Figure  3 is  a limited  set  for  the  carbon/hydrogen/oxygen  species.  All 
reactions  are  not  shown  in  these  diagrams,  and  reactions  linking  the 
families  are  not  all  included. 


COMPUTATIONAL  RESULTS 

The  calculations,  which  produced  the  selected  sample  results  presented 
below,  were  fully  time-dependent  and  diurnal.  Many  calculations  of  this 
type  use  a number  of  equilibrium  relationships  and/or  a steady-state 
solar  input,  and  are  therefore  less  realistic.  Although  the  calculations 
were  diurnal,  and  the  results  obtained  are  available  for  any  time  (day, 
night,  or  transition),  only  the  results  for  noontime  are  presented.  The 
following  figures  are  therefore  the  model -derived  characteristics  of  the 
stratosphere  at  noon,  in  late  September  at  latitudes  near  32°  N. 

The  categories  of  calculated  results  chosen  for  presentation  include: 

(1)  particle  densities  of  selected  important  species,  (2)  mixing  ratios 
of  these  constituents,  (3)  rates  of  several  sets  of  important  photolytic 
and  chemical  reactions,  (4)  total  formation  and  removal  rates  of  some 
key  species,  and  (5)  equivalent  transport  contribution  to  a few  molecular 
densities.  All  are  presented  as  a function  of  altitude  between  10  and  50 
kilometers.  Other  types  of  information  can  be  derived  from  the  model, 
and  equivalent  results  for  other  latitudes  and  periods  of  the  year  are 
readily  calculable. 


COMPOSITION-PARTICLE  DENSITIES 

Figures  4a  to  4k  are  the  calculated  altitude  profiles  of  the  particle 
densities  of  ground-state  oxygen  atoms  0 (^P),  excited  oxygen  molecules 
02(^4),  ozone  (O3),  nitric  oxide  (NO),  nitrogen  dioxide  (NO2),  nitric 
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acid  (HNO3),  nitrous  oxide  (N2O),  water  (H2O),  carbon  monoxide  (CO), 
hydroxyl  radical  (OH)  and  methane  (CH,,).  (Units  are  cm“^.)  All  of  these 
species  ere  of  fundamental  importance  in  stratospheric  chemical  kinetics, 
and  all  have  been  measured  in  the  stratosphere.  However,  the  number  of 
measurements  is  very  limited  except  for  ozone  and  water,  and  even  with 
these  two,  variability  and  uncertainty  are  high. 

Figures  4f  to  4p  are  calculated  particle  densities  of  some  of  the 
speculative  species  which  have  not  yet  been  detected  in  the  stratosphere, 
but  whose  presence  is  required  by  photolytic  and  chemical  considerations. 
These  species  include  electronically  excited  oxygen  atoms  0(^0),  perhy- 
droxyl  radicals  (HO;),  hydrogen  peroxide  (H2O2),  formaldehyde  (HCHO)  and 
nitrous  acid  (HNO2).  In  the  absence  of  measurements,  the  only  informa- 
tion available  on  such  constituents  comes  from  models  and  related  theory. 
In  particular,  the  0(^D)  atoms  and  the  HO2  radicals  are  involved  in  some 
critically  important  atmospheric  reactions,  including: 


0(^D) 

+ H2O  OH  + OH 

(73)* 

0(^D) 

+ N2O  . NO  + NO 

(71) 

HO2  + 

OH  > O2  + H2O 

(132) 

COMPOSITION-MIXING  RATIOS 

Figures  5a  to  5p  present  the  same  data  as  4a  to  4p  expressed  as  mixing 
ratios:  the  particle  density  of  the  species  divided  by  the  total 

particle  density  of  all  the  species  present  in  the  volume  considered. 
This  ratio  of  numbers  of  particles  is  termed  the  volume  mixing  ratio. 

The  ratio  of  the  masses  (mass  mixing  ratio)  is  also  frequently  used, 
particularly  for  water.  In  these  figures,  the  abscissa  designation  10“® 
is  one  part  per  million  (1  ppmv),  and  10"'^  is  one  part  per  billion 
(1  ppbv). 

The  figures  provide  no  new  basic  information,  but  allow  direct  consid- 
eration of  the  fraction  of  the  atmosphere  that  a given  constituent 
represents.  They  emphasize  the  fact  that  species  which  play  dominant 
roles  in  many  atmospheric  processes  are  present  in  only  extremely  small 
proportions . 


RATES  OF  CHEMICAL  AND  PHOTODISSOCIATION  REACTIONS 

To  understand  atmospheric  characteristics  and  behavior,  it  is  necessary 
to  know  what  processes  are  occurring  and  the  role  and  relative  impor- 
tance of  each  under  the  range  of  conditions  encountered.  As  indicated 
previously,  the  processes  are: 


*The  reaction  numbers  used  throughout  are  as  used  in  [2]. 
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Radi ati ve  - concerned  with  the  transmission,  scattering,  absorption 
and  emission  of  energy  in  the  infrared,  visible,  and  ultraviolet  regions, 
and  the  thermodynamic  consequences  on  atomic/molecular  scales, 

Dynami c - covering  the  range  of  particle  motions  from  molecular  to 
glob^  scales . 

Chemical  Kinetic  - treating  the  multitude  of  chemical  reactions 
among  the  many  species  present. 

A useful  method  of  deriving  and  demonstrating  the  role  of  various  chemical 
reactions  is  to  calculate  the  rate  at  which  a given  reaction  proceeds 
under  the  conditions  considered.  Figures  6a  to  6n  are  the  calculated 
altitude  dependencies  of  the  rates  of  a number  of  important  reactions 
for  the  noon,32'"N,  September  situation.  The  reactions  are  listed  in 
Table  2.  Some  of  the  reactions  are  grouped  to  allow  direct  comparison  of 
certain  critical  relationships.  Since  the  relevance  or  role  of  a specific 
reaction  or  set  of  reactions  is  dependent  on  the  particular  problem  under 
consideration,  the  following  discussion  will  not  be  in  the  context  of 
stratospheric  pollution,  weapons  effects,  meteorology,  or  other  specific 
present  concerns,  but  will  be  limited  to  a few  brief  general  remarks. 

Figure  6a  shows  the  rates  for  the  most  important  atmospheric  reaction: 
the  initial  production  of  oxygen  atoms  by  photodissociation  of  oxygen 
molecules,  which  sets  the  stage  for  most  of  the  stratospheric  chemistry 
of  consequence.  It  should  be  noted  how  rapidly  this  photolysis  falls 
off  below  25  kilometers  altitude  - the  rate  at  10  kilometers  is  less 
than  one  thousandth  that  at  25  kilometers. 

Figure  6b  represents  the  two  principal  modes  of  photodissociation  of 
ozone.  In  reaction  5 the  products  are  ground-state  atoms  and  molecules, 
while  reaction  7 produces  electronically  excited  (and  more  reactive) 
atoms  0(^D)  and  molecules  02(^a).  Reaction  5 dominates  below  30  kilo- 
meters, and  reaction  7 above.  Ozone  formation  (reactions  36  + 37)  is 
shown  in  Fig.  6c;  and  the  comparison  of  this  three-body  formation 
rate  with  the  photolytic  destruction  (the  sum  of  reactions  5 and  7)  is 
presented  in  Fig.  6d.  It  should  be  noted  that  these  reaction  rates 
are  very  large  compared  with  other  atmospheric  reactions  and  that  they 
are  almost  equal  at  all  altitudes,  so  that  an  equilibrium  among  0,  O2 
and  O3  closely  follows  the  solar  flux  input. 

Three  important  ways  in  which  the  oxygen  atoms  (produced  in  Fig.  6a  and 
brought  somewhat  into  equilibrium  in  Fig.  6d)  enter  into  the  chemistry 
is  demonstrated  in  Fig.  6e.  These  are  three  of  the  principal  removal 
mechanisms  of  oxygen  atoms.  The  removal  of  0 and  O3  by  the  much- 
publicized  NO^  "catalytic"  cycle  is  shown  in  Fig.  6f.  The  O+NO2  reaction 
(43)  is  much  slower  over  this  altitude  range  and  is  therefore  the  rate- 
controlling reaction  for  this  "ozone  depletion"  process. 
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''  TABLE  2 

REACTIONS  SELECTED  TO  ILLUSTRATE  ALTITUDE-DEPENDENCE 
OF  RATES  FOR  FIGURES  OF  SECTION  3 


(1) 

O2  + hv 

-V 

0 + 0 

(5) 

O3  + hv) 

- 

0 + O2 

(7) 

O3  + hv 

- 

0(1D)  + 02(1a) 

(11) 

NO2  + hv 

0 + NO 

(32) 

HNO 3 + hv 

'V 

OH  + NO2 

(36) 

0 + O2  + O2 

- 

O3  + O2 

(37) 

0 + O2  + N2 

-V 

O3  + N2 

(38) 

0 + O3 

->- 

O2  + O2 

(43) 

0 + NO2 

O2  + NO 

(50) 

0 + HO2 

->• 

C2  + OH 

(71) 

0(1D)  + N2O 

->■ 

NO  + NO 

(73) 

0(1D)  + H2O 

OH  + OH 

(94) 

O3  + NO 

•y 

O2  + NO2 

(95) 

O3  + NO 

-y 

O2  + NO2* 

(101) 

O3  + OH 

•> 

O2  + HO2 

(115) 

NO  + NO3 

-> 

NO2  + NO2 

(117) 

NO  + HO2 

~y 

OH  + NO2 

(119) 

NO  + HO2  + M 

- 

HNO3  + M 

(123) 

OH  + NO2  + M 

y 

HNO3  + M 

(126) 

H + O2  + M 

■> 

HO2  + M 

(132) 

OH  + HO2 

•V 

O2  + H2O 

(134) 

OH  + CO 

>■ 

H + CO2 

(135) 

OH  + CH4 

■V 

H2O  + CH3 
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Also  important  in  the  NO/0;  problem,  and  the  utlimate  controlling  factor 
when  combined  with  HNO3  dynamic  removal,  is  the  three-body  formation  and 
photolytic  loss  of  nitric  acid  in  reactions  123  and  32  as  plotted  in 
Fig.  6g.  The  balance  among  formation,  loss,  and  transport  removal  of 
HNO3  will  help  establish  the  NO2  density  since  it  will  remove  NOj  from 
the  equilibrium  concentration  that  the  reactions  in  Fig.  6h  tend  to 
establish.  These  are  the  rates  for  formation  of  NO2  by  the  NO+O3  reaction 
and  for  its  photodissociati ve  destruction.  Below  35  kilometers,  these 
are  very  nearly  equal. 

Figure  6i  compares  trie  NO+O3  rate  with  other  important  nitric  oxide 
reactions,  showing  its  dominance  at  all  altitudes.  Figure  6j  demon- 
strates the  role  of  the  excited  oxygen  atoms  0(^D)  produced  by  ozone 
photolysis  as  the  principal  primary  source  of  OH  radicals  and  NO  molecules. 

Figure  6k  gives  the  rates  for  three  of  the  most  important  reactions  of 
these  hydroxyl  radicals,  while  Fig.  6i  shows  their  role  in  destruction 
of  the  natural  and  anthropogenic  carbon  monoxide  (CO)  and  methane  (CH..) 
that  move  into  the  stratosphere  from  the  earth's  surface. 

Figure  6m  shows  the  principal  formation  mechanisms  of  the  perhydroxyl 
radical  {HO2).  Figure  6n  compares  formation  of  nitric  acid  from  NO  and 
from  NO2  showing  the  potential  importance  of  the  OH/NO2  process  in  the 
lower  stratosphere. 


TOTAL  FORMATION  AND  REMOVAL  RATES 

Figures  7a  and  7b  show  the  calculated  noontime , 32° N latitude,  September 
altitude  dependencies  of  the  formation  and  removal  rates,  by  all  appli- 
cable chemical  and  photolytic  processes,  for  eight  minor  constituents: 

0(3P),  0(iD),  O3,  NO,  NO2,  OH,  HO2  and  HNO3 

The  formation  rate  is  essentially  equal  to  the  removal  rate  for  all  but 
HNO3  in  this  near-equilibrium  noontime  situation.  Hence,  for  all  but 
HNO3  only  one  curve  is  seen  for  both  formation  and  removal.  This  was 
implied  for  ozone  in  Fig.  6d  since  the  entire  production  of  ozone  is  by 
reactions  36  and  37,  and  almost  all  of  the  loss  is  by  the  photodissociation 
reactions  5 and  7. 

Figure  7a  also  demonstrates  that  oxygen  atom  formation-removal  rates  are 
virtually  equal  to  those  of  ozone  in  the  altitude  range  considered. 
Similarly  the  NO  and  NO2  pair  have  essentially  equal  production-loss 
rates.  The  hydroxyl  radical  (OH)  format  ion- removal  rates  shown  in 
Fig.  7b  are  approximately  equal  to  those  of  the  perhydroxyl  radical 
(HO2)  in  the  lower  stratosphere,  but  are  a small  factor  larger  in  the 
upper  stratosphere. 
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TRANSPORT  CONTRIBUTION  TO  PARTICLE  DENSITIES 

The  continuity  equations  used  to  describe  the  time  dependencies  of  the 
particle  densities  (n^. ) of  the  atmospheric  chemical  constituents  contain 
chemical  formation  terms  (F^),  chemical  and  photolytic  loss  terms  (L^-)> 
and  a divergence  term  to  describe  the  net  gain  or  loss  of  the  particles 
due  to  transport  through  the  volume  considered: 

^ [n.]  = F.  - L.  - vpv 

The  effect  of  the  great  variety  of  transport  processes  on  atmospheric 
composition  is  very  poorly  understood,  and  is  usually  incorporated  into 
models  in  somewhat  artificial  fashions.  One  purpose  of  the  present  model 
is  to  reverse  this  procedure,  since  the  photolytic  and  chemical  processes 
have  been  much  better  characterized,  and  derive  mean  characteristics  of 
the  vertical  contribution  to  transport  using  atmospheric  measurements 
combined  with  detailed  chemical  kinetic  treatment.  A prel imi nary  part 
of  this  study  is  presented  partially  and  briefly  in  Figs.  8a  to  8e, 

These  figures  show  the  regions  of  the  stratosphere  where  the  net  result 
of  the  transport  through  the  volume  is  an  increase  of  the  particle  den- 
sity (points  plotted  +)  or  a decrease  (points  plotted  -).  Alternatively, 
this  can  be  described  as  a region  where  the  chemical/photolytic  loss 
rate  is  greater  than  the  chemical  formation  rate  (+),  or  where  the 
formation  rate  is  greater  (-).  The  constituents  considered  are  0i(,  NO25 
N2O,  H2O,  and  HNO3. 

In  Fig.  8a,  in  the  region  above  30  kilometers  altitude,  removal  of  ozone 
by  photolysis  and  reaction  with  NO,  0,  OH,  etc.,  is  faster  than  the  three- 
body  chemical  formation,  and  ozone  is  transported  into  the  region.  Below 
the  transition  at  30  kilometers,  there  is  a net  formation  and  ozone 
transport  is  out  from  the  region.  The  units  of  net  loss  or  gain  are 
cm“^  sec"i. 

For  NO2  (Fig.  8b),  excess  formation  is  above  17-19  kilometers,  with 
deficient  formation  below  and  a net  flow  into  the  lower  region.  Nitrous 
oxide  (N2O)  has  a net  chemical/photolytic  loss  throughout  the  strato- 
sphere, and  the  equivalent  contribution  to  the  N2O  density  by  the  upward 
flow  is  at  the  rates  shown  in  Fig.  8c. 

The  reverse  of  this  occurs  for  water  (Fig.  8d)  with  chemical  formation 
exceeding  losses  through  the  10-50  kilometer  range.  Under  the  conditions 
considered,  nitric  acid  losses  dominate  above  20  kilometers,  with  a flow 
upward  and  downward  from  the  10-20  kilometer  region  where  there  is 
excess  formation,  as  indicated  in  Fig.  8e. 

The  results  presented  in  Figs.  8a  through  8e  are  exploratory  and  pre- 
liminary. 


Results  of  some  of  the  experimental  measurements  on  STRATCOM  VI  reported 
in  [3]  have  been  plotted  on  the  graphs  of  calculated  particle  densities 
and  mixing  ratios.  These  include  ozone  (Figs.  4c  and  5c),  nitrous  oxide 
(Figs.  4g  and  5g),  carbon  monoxide  (Figs.  4j  and  5j)  and  methane 
(Figs.  4k  and  5k).  It  should  be  noted  that  since  the  model  uses  mean 
values  as  inputs,  the  results  are  plotted  as  smooth  curves.  The  actual 
stratosphere  can  have  sharp  variations  of  a parameter  over  a small 
altitude  interval,  and  an  equivalent  plot  of  measured  altitude  dependence 
could  be  a jagged  line. 

> 

*'  However,  even  with  this  consideration,  the  relatively  good  agreement  of 

’ the  computational  results  with  the  measurements  for  the  given  season, 

time,  latitude,  and  altitude,  gives  some  degree  of  confidence  in  the 
validity  of  the  model  and,  consequently,  in  the  other  results  calculated. 

j Similarly,  some  earlier  results  of  J.  G.  Anderson,  at  the  same  latitude, 

for  the  short-lived  highly  reactive  oxygen  atoms  and  hydroxyl  radicals 
are  plotted  comparatively  in  Figs.  4a  and  4i , respectively.  The  oxygen 
atom  results  are  for  November  1974,  10:30  a.m.  [6].  The  length  of  the 
line  indicates  the  experimental  uncertainty.  For  the  hydroxyl  radical 
f (OH),  the  dotted  line  is  for  measurements  in  July  and  the  solid  lines 

i for  January  1976  [7]. 


Figure  1.  Diagrammatic  representation  of  reactions  among  atmospheric  constituents  composed 
only  of  oxygen  and  nitrogen  atoms. 


Fi|ur«  4.  Composition-particlt  densitiss. 


N,0  MIXING  RATIO  H,0  MIXING  RATIO 


(lOl)  t OH  O2  + HO2 
023;  M02+-OH-!^  HNO3 


(u'<)  Haniiiit 


HNO,  f'LUX  DlV[RGtNCr  (CM'-^SfC 
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